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ABSTRACT 

It has been suggested that the high metallicity generally observed in active galactic nuclei (AGNs) and 
quasars originates from ongoing star formation in the self-gravitating part of accretion disks around the su- 
permassive black holes. We designate this region as the star forming (SF) disk, in which metals are produced 
from supernova explosions (SNexp) while at the same time inflows are driven by SNexp-excited turbulent vis- 
cosity to accrete onto the SMBHs. In this paper, an equation of metallicity governed by SNexp and radial 
advection is established to describe the metal distribution and evolution in the SF disk. We find that the metal 
abundance is enriched at different rates at different positions in the disk, and that a metallicity gradient is set up 
that evolves for steady-state AGNs. Metallicity as an integrated physical parameter can be used as a probe of 
the SF disk age during one episode of SMBH activity. In the SF disk, evaporation of molecular clouds heated 
by SNexp blast waves unavoidably forms hot gas. This heating is eventually balanced by the cooling of the hot 
gas, but we show that the hot gas will escape from the SF disk before being cooled, and diffuse into the BLRs 
forming with a typical rate of ^ IMoyr"'. The diffusion of hot gas from a SF disk depends on ongoing star 
formation, leading to the metallicity gradients in BLR observed in AGNs. We discuss this and other observable 
consequences of this scenario. 

Subject headings: black hole physics — galaxies; evolution — quasars: general 



1. INTRODUCTION 

It is widely accepted that active galactic nuclei (AGNs) and 
quasars are powered by accretion onto supermassive black 
holes (SMBHs) (Lynden-Bell 1969). The standard accre- 
tion disk model (Shakura & Sunyaev 1973) or its variants 
that account for the so-called "big blue bump" (Shields 1978; 
Malkan 1983; Laor & Netzer 1989; Sun & Malkan 1989; 
Brunner et al. 1997; Brocksopp et al. 2006; Lu 2008) show 
that a typical accretion rate is a fraction of the Eddington limit. 
Over the last four decades, much attention has been paid to the 
radiation process and dynamics of the accretion disk, but far 
less is known about the mechanism driving such an accretion 
rate. There is little understanding of how the SMBHs are fed 
at the ^ 1 pc scale, or about the outer boundary condition of 
the Shakura-Sunyaev accretion disk. 

Two important clues can help resolve these wide-open is- 
sues. The first is the role of the vertical self-gravity of the 
outer part of the accretion disk. Paczynski (1978) pointed out 
that quasar accretion disks are so massive that self-gravity 
dominates. As a consequence, the disk unavoidably col- 
lapses into clumps, among which collisions provide a natu- 
ral mechanism for viscosity. This was recently stressed again 
by Rafikov (2009). However, Kolykhalov & Sunyaev (1980) 
suggested that star formation cannot be avoided in such a mas- 
sive disk. Since this pioneering work, many attempts have 
been made to construct models that incorporate star formation 
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in the disk as well as take account of its roles in the viscosity 
and related issues (Sakimoto & Coroniti 1981; Lin & Pringle 
1987; Shlosman & Begelman 1987; Shlosman et al. 1989; 
Shlosman & Begelman 1989; Gammie 2001; Goodman 2003; 
ColUn & Zahn 1999, 2008; Thompson et al. 2005; Tan & 
Blackmann 2005; Nayakshin & Sunyaev 2005; Nayakshin et 
al. 2008; and the recent numerical simulations of Hobbs et al. 
2011; Jiang & Goodman 2011). All of these studies reach a 
similar conclusion: star formation in the dense accretion disk 
is so efficient that the remaining gas cannot adequately feed 
the SMBHs unless there is an external torque to transport an- 
gular momentum (e.g. Thompson et al. 2005). 

The second major clue is the high metal abundances ob- 
served in quasar broad-line regions (BLRs), showing that 
there is a fast process of metal enrichment in the nuclear re- 
gions of AGNs unless the enrichment has been done before 
reaching the nuclear regions (Hamann & Ferland 1993; see a 
review of Hamann & Ferland 1999; Collin 1999). Generally, 
the metallicity of host galaxies is significantly lower than that 
of the BLR, with even the circumnuclear regions of galax- 
ies having only slightly more than solar metallicity (Dors et 
al. 2008). The BLR metallicity con-elates with SMBH mass 
or AGN luminosity (Hamann & Ferland 1999; Warner et al. 
2003; Matsuoka et al. 2011). Recently, Simon & Hamann 
(2010) found that the BLR metallicity does not correlate with 
star formation rates of host galaxies as traced by far-infrared 
emission. These observational indications provide evidence 
for intense star formation in a massive disk associated with 
the central engine rather than that the BLR metals are trans- 
ported from the host galaxies. However, the problem remains 
that SMBHs must be fed with a sufficiently high mass rate 
to produce the observed luminosity. It has seemed irrecon- 
cilable for AGN to simultaneously have both high metallicity 
and large accretion rates. 

A solution to this puzzle is to employ turbulence excited by 
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supernovae explosions (SNexp) of massive stars during star 
formation as a means to transport angular momentum (Chen 
et al. 2009; Wang et al. 2010; Hobbs et al. 2011). Recently, 
Wang et al. (2010) stressed the roles of star formation in metal 
production in the inflows developed from the inner edge of the 
dusty torus, driving not only a sub-Eddington-limit accretion 
rate of the SMBHs, but also the formation of a compact nu- 
clear star cluster after multiple AGN episodes. They find that 
SNexp have two main effects in the SF disk: 1) they deliver 
gas inward, reducing local star formation rates; and 2) the top- 
heavy initial mass function quickly feeds processed gas back 
into the gaseous disk, enhancing the metallicity. The gas sup- 
ply to the SMBHs is thus maintained with a rate typical of 
quasars, while at the same time the metal-rich properties such 
as the metallicity-luminosity trend can be quantitatively ex- 
plained (Wang et al. 2010). The BLR metallicity measured 
from the N V /C IV ratio (Hamann & Ferland 1993) directly 
reflects the metallicity in the SF disk as long as the BLR con- 
tinues to be fed from the SF disk. Quantitative predictions of 
the properties of the SF disk and the way in which it feeds 
the BLR are urgently needed so that observational tests can 
be made in the future. 

This paper is structured as follows. §2 gives the model of 
SF disk in light of Toomre's parameter and justifications for 
it. §3 is devoted to building a metallicity equation for the SF 
disk that can predict its metallicity gradient and evolution. §4 
focuses on the development of thermal diffusion of hot gas 
from the SF disk. We show that the diffusion rate can reach as 
high as a fewMoyr"' in quasars. Observational consequences 
are then discussed in §5. Conclusions and some further brief 
discussion are given in the last section. 

2. THE GENERAL STRUCTURE OF THE STAR FORMING DISK 

2.1. Underlying physical processes in the disk 

We start by summarizing the results published to date about 
star forming disks. Most of the release of gravitational energy 
around an SMBH occurs in a region within about 3 - 100/?sch 
of the black hole, where /?sch = 2GM,/c^ is the Schwarzschild 
radius, G is the gravity constant, M, is the SMBH mass and c 
is the light speed. This part of the accretion disk is understood 
much better than its outskirts. To minimize the uncertainties 
due to the opacity etc., we start with the Shakura-Sunyaev 
disk for the self-gravity radius given by equation ( 1 8) in Laor 
&Netzer(1989) 



R.c,^ 511 a^l'M-^"m'l'>R, 



Sch, 



(1) 



where a is the viscosity parameter, rh = M/MEdd, A^Edd = 
2.22?7o\M8M0yr"', Mg =M,/1O**M0 and ryo.i = ?7/0.1 is the 
radiative efficiency, assuming Toomre's parameter 2=1. This 

is very similar to /?sg « 652 a'*/^' fi-^^^^J'' e^J^^ Alf'''^' Rs,b, 
where kr is the Rossland opacity, /i is the mean molecular 
weight, i?Edd = ^Boi/^Edd IS the Eddington ratio and LboI is the 
bolometric luminosity (Collin & Hure 1999). Beyond this ra- 
dius, the disk is unstable to self-gravity, inducing star forma- 
tion. In the self-gravitating regions, the radiation from star 
formation is given by L, = rj/IZc^, where TZ is the star for- 
mation rate and 77* is the efficiency of converting mass into 
radiation through star formation. In the same regions, the dis- 
sipated gravitational energy bound by the SMBH is given by 
^acc ~ ^'^Sch^^c^. Comparing the two energy sources, we 
find that stellar energy dominates over the gravitational at a 

critical radius i?c ~ ri*(M /Ti)Rsch = 10-'??-3(M/7t)i?sch, where 



?7_3 = ?7,/10 -^ (?7* = 10 ^ for a Salpeter initial mass function 
in Thompson et al. 2005). This clearly states that Rsc ^ Re 
for a SF disk with TZ ^ M, implying that the self-gravitating 
regions are governed by stellar energy rather than the gravity 
of accreting gas bound by SBMHs. 

The underlying physics mainly involves several dynami- 
cal processes coupling with thermal ones. Fragmentation 
into clouds happens when the cooling timescale fcooi < 3il~', 
where ft is the rotation frequency (Gammie 2001). For a Kep- 

lerian rotating disk, O"' ?» 15.0 r/ yr, where r4 =/?/10^/?sch, 
and fcooi = n„2 kT/Au^in^^ , T) w 5.5 T^^'^yr, where A(n„^ , T) = 
4.2 X 10~-"n„ r-'-'erg s"' cm"^ is the cooling function and hh, 
is the density of hydrogen molecules (Smith & Mac Low 
1997). Here estimates of fcooi are conservative since we do 
not include the metallicity. We find that the fragmentation of 
the disk is inevitable, leading to star formation. 

The self-gravitational instability drives collisions among 
clouds, leading to effective viscosity (e.g. Duschl & Britsch 
2006), but it could be less important than the viscosity 
due to turbulence produced by the SNexp. The disper- 
sion velocity obtained from the fragmentation is given by 
y^i^p = (GMdO)i/3(f<;oii/fgrav)'''^ (ColHn & Zahn 2008, here- 
after CZ08), where M^ is the cloud mass, and fcoii and fgrav are 
the collision and gravitational interaction between two clouds, 
respectively. For massive clouds with a mass of IO-'Mq, we 
have Vdisp « 156 r^ km s"', which is too small compared 

with the Keplerian rotation VKep = 2.1 x lO-' r^ ' km s"' to 
excite strong enough turbulence for transportation of angular 
momentum in the disk (< 1 pc). The SNexp-driven turbulent 
velocity is of the order of the Keplerian velocity (see §3.1) and 
is regarded as the main sources of turbulence, that is to say, 
q; « 1 is suggested by 2D numerical simulation (Rozyczka 
et al. 1995). We note that CZ08 underestimate the SNexp- 
driven viscosity since the SNexp would be anisotropic in the 
fast rotating disk (see detailed arguments in §3.1). 

Detailed processes for forming stars in the disk were stud- 
ied by ColHn & Zahn (1999, hereafter CZ99) and CZ08. They 
considered the following processes: accretion from the SF 
disk to form the stars, stellar winds subsequently shutting 
off such accretion flows, the effects of the accretion process 
opening gaps in the SF disk around the forming stars, and 
migration of stars toward the SMBH. CZ99 argue that these 
stars may be driven to migrate toward the black hole by den- 
sity waves in the disk, but may be influenced by the cavity 
opened by the stellar wind. If the timescale of migration to- 
ward the SMBH is shorter than the evolutionary timescale of 
main sequence stars, the SNexp-driven turbulence does not 
work. However, there is actually another mechanism to stop 
migration toward the SMBHs, which is that the accretion of 
stars from the disk increases not only its mass, but also the 
angular momentum of the stars. The gain of angular momen- 
tum makes the stars migrate outward. The accretion timescale 
of stars is face ~ lO'* (m*/lOM0)"'yr inside the disk environ- 
ment (Artymowicz et al. 1993). Comparing face with tj^igr 
given by equations (28) and (29) in CZ99, we find that the 
stars stay at approximately the location where they were born 
since generally face < fmigr works, where f^igr is the radial mi- 
gration timescale. This allows us to neglect the migration of 
stars, without obvious influence on the evolution of the stars. 

Finally, stars reach their maximum mass when stellar winds 
balance their accretion from the disk (CZ08) and they then 



Star Formation in Self-gravity Accretion Disks in AGNs 



maintain that mass until they explode as SN. The complicated 
interplay between these processes leads to uncertainties in the 
maximum and minimum masses of stars. The initial mass 
function of stars is not obtained in a self-consistent way, but 
the uncertainties arising from accretion and stellar winds re- 
main as free parameters in the initial mass function (see §3.2). 



man (2003) 




R (pc) 



Fig. 1. — Structure of the self-gravitating disk, from Eqs (2), (3) and (4). 
Here we use M. = W^Mq , k = Kes = 0.4cm^g"' , a = 0. 1 and ^Edd = 0.3. The 
surface density rate of star formation is very similar to the numerical results 
given by Wang et al. (2010). For simplicity, dependence of the opacity on the 
metallicity is neglected in the calculations. 



2.2. The model of a star forming disk 

Wang et al. (2010) computed numerical models of the 
structures of star forming disks. In order to carry out an ana- 
lytical study of the effects of the SF disk, we use a simplified 
version of its structure, which are approximately the simple 
version of a self-gravitating disk (e.g. Goodman 2003). The 
smoothed structure of the self-gravitating disk is determined 
by Toomre's parameter Q, and is justified in this section. It 
has been demonstrated by Thompson et al. (2005) that such a 
disk is supported by the radiation pressure due to star forma- 
tion, yielding Q=\. The surface density Eg^s and the height of 
the SF disk Ha are given by equations (18) and (20) in Good- 



Egas= 1.6x10' 



tEdd 



OiQAm.l 



1/3 



m; 



-2/3 -3/2 



'3 



Mopc" 



and 



//d = 4.6x10"^ 



^Edd 
"0.1 '70.1 



1/3 



Mt'WI^ 



pc, 



(2) 



(3) 



where r-}, = /?/ 10^/?sch- The validity of this smooth model of a 
star forming disk is discussed below. 

Fragmentation of the disk depends on the rotation, den- 
sity, temperature, metallicity and cooling of individual clouds, 
whose necessary condition is generally fulfilled as discussed 
in the previous section. Details of the star formation process 
and macrophysics are poorly understood (McKee & Ostriker 
2007), especially in the accretion disks. Here we will use 
the empirical Kennicutt-Schmidt law of star formation rates 
(Kennicutt 1998). In NGC 5194 (Kennicutt et al. 2007) 
and NGC 4252 (Rahman et al. 2010), this law works over 
all scales from the entire galaxy to the nuclear region, al- 
though in M33 (Onodera et al. 2010) it does appear to break 
down at the scales of giant molecular clouds. As a way to 
move ahead and test the observational consequences of star 
formation, we will assume that the Kennicutt-Schmidt law 
is at least approximately correct for the accretion disk. The 
results found here will still be qualitatively correct unless 
the star formation rates are completely unrelated to the gas 
surface density. According to the Kennicutt-Schmidt law, 
S* = 2.5 X 10"^Sg;,|MQyr"'kpc"^, where Sgas is the surface 

density of gas in unit of Mq pc~^, and we have the surface 
density rate of star formation 



S* = 1.9x10' 



^F 



ao.1'70.1 



0.47 



m; 



-0.93-2.1 



rj^' Meyr-'pc-^ (4) 



We find that S* depends more strongly on the SMBH mass (as 
M;°^') and the distance to the SMBHs than on the Eddington 
ratio (as shown in Figure 1). This leads to the conclusions that 

1) AGN with less massive SMBHs may have different proper- 
ties from ones with more massive SMBHs, yielding observa- 
tional implications such as narrow line Seyfert 1 galaxies; and 

2) there will be a strong S* gradient within individual objects 
that will lead to a steep metallicity gradient in the SF disk 
though there is radial advection as shown in §3.1, resulting in 
a difference in the metallicity measured by different indica- 
tors. We discuss these two points in the following sections. 

The clumpy SF region is composed of molecular clouds, 
evolving stars, SNexp (including the SN shells), and diffuse 
hot gas. The structure of the disk is determined by the mu- 
tual interaction among these components (McKee & Cowie 
1975). For simplicity, we use a smooth model of a self- 
gravitating disk to represent in analytical form the averaged 
properties of the disks, but we explicitly consider a clumpy 
structure when we determine the way in which gas is heated 
by SNexp to diffuse into the broad line regions. The va- 
lidity of the "smooth" thin disk can be justified by the fact 
that the characteristic distance (da) among clouds in the star 
forming regions is much less than the global size of the disk. 

The self-gravitating timescale fcoUapse = {Gpc) , where pc is 
the density of clumps, should not be longer than the crossing 
timescale of turbulence through a cloud ftur = 'c/Vtur, where 
Ic and Vtur are the characteristic length of clouds and turbu- 
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lence velocity. Otherwise, the cloud will be destroyed to form 
stars. This yields the density of clouds as pc = G~' (Vtui/^e) ■ 

We then have d^\ = (pc/p) h- = ^tui- (o^H/Gp) , where p 
is the mean density of the disk and l^ « aH is used. For typ- 
ical values of the disk, p = 10~'^g cm"-', Vtuj = lO-'km s~' and 
H = lO'^cm at 7? = lO'^TJsch, and we have d^i « 2.5 x lO'^cm, 
which is significantly smaller than the radius of the disk. The 
"smoothed" model is a quite good approximation of the star 
forming dislfl 

We should point out that a constant opacity is used for the 
self-gravitating disk as by Goodman (2003). Collin & Hure 
(1999) study the influence on the disk of opacity, as a func- 
tion of density and temperature. They show that the surface 

density of the disk is given by Sgas ex Kj/ and Sgas oc k,/ for 
solar abundance and zero-metallicity, respectively (see equa- 
tions 8 and 12 in Collin & Hure 1999). Though the opacity 
kr can change by a factor of 10^, approximating the opacity 
as a constant only leads to a factor of a few uncertainty in the 
surface density. We note that the disk structure displays sev- 
eral solutions with thermal instability etc. depending on the 
opacity (Hure et al. 1994). The resulting star formation in the 
disk is not very sensitive to the opacity based on the simplified 
model of Goodman (2003), but we keep these approximations 
in mind. 

Finally, we note that star formation in the disk is a self- 
regulating process in light of the Toomre's parameter, as ar- 
gued by CZ99. In this paper, we neglect this dependence of 
star formation, but use the Kennicutt-Schmidt law for simplic- 
ity. Although the Kennicutt-Schmidt law might be modified 
somehow in the disk, this approach still illustrates the effects 
of star formation there. We would like to point out that star 
formation takes place in the regions with a Keplerian veloc- 
ity of (1 -3) X lO^km s"^ Detailed process of star formation 
needs an efficient way to get rid of the angular momentum of 
the protostars, perhaps, magnetic fields as a potential mecha- 
nism could be invoked to do so. 

3. METALLICITY WITHIN THE STAR FORMING DISK 

3.1. SNexp evolution in the SF disk 

In the SF disk, the evolution of a SNexp (here meaning 
the SN shell) fully depends on its surroundings. It may con- 
veniently be divided into three stages (e.g. Wheeler et al. 
1980 and Shull 1980): 1) rapid ejection followed by expan- 
sion with a nearly constant velocity; 2) adiabatic expansion in 
the Sedov-phase; 3) the snowplow phase governed by radia- 
tion loss when a dense shell forms at the outer boundary. A 
SNexp in the SF disk is different from the well-known case 
for a homogeneous medium, because the SF disk is much 
denser so that "blowouts" in the direction perpendicular to 
the disk will be very rare. Rozyczka et al. (1995) studied 
the evolution of SNexp in a disk within a Ipc scale through 
2-D numerical simulations, and found that the SNexp can in 
principle provide turbulent viscosity of a ~ 0.1. Here we 
only give a rough estimate of the SNexp evolution in the SF 
disk. According to Shull (1980), for an isotropic SNexp, its 
Stage 1 ends when the swept mass is comparable with the 
ejecta of the SNexp, yielding the radius of the SNexp at this 

phase 7?si = 4.1 x 10"^n^4 M .Qpc, where n\^ = «/10''*cm"-' 

and Mej,o = M^jKf'MQ is the ejected mass. The duration of 

' We note that the SF disk could be influenced by the stars for the extreme 
case that the mass of stars is very large. 



-■X -1/3 

Stage 1 is given by fsi = 4.0 x 10 ni4 yr if the initial ejec- 
tion velocity is lO^km s"'. It is found that Rsi < //d and 

hi ^ fK, where fK = 27r/51 = 8.8 Mg'rj yr is the timescale 
of the Keplerian rotation. Stage 1 of a SNexp is not able 
to be break out of the SF disk (Rozyczka et al. 1995). 
Wheeler et al. (1980) demonstrate that the Sedov-phase does 
not exist provided the density of surrounding gas is larger 

than 5 x IQ^El^l (Mej/IOM©) cm"^, which is the case here. 
Clearly, the SNexp directly enters Stage 3 without radiation 
loss. Most of the kinetic energy of the SNexp will be chan- 
neled into the kinetic energy of clumps. The maximum ra- 
dius of the SNexp is typically R^n ^ Rs\ ^ 0.1//d, implying 
that the viscosity driven by the SNexp will be a '-^ 0.1 - 1, 
which is consistent with numerical simulation of Rozyczka et 
al. (1995). 

On the other hand, the kinetic energy of the SNexp heats the 
surface of clumps through blast shocks when the SNexp hap- 
pens. As we show next, the SNexp ensures that evaporation of 
molecular clouds via shocks can drive thermal diffusion of the 
hot gas above the SF disk into the broad line regions. Table 
1 gives a summary of the present model including the typical 
parameter values for an SMBH with IO'^Mq, which are useful 
for future observational tests. 

There are debates as to the strength of SNexp-driven viscos- 
ity in the star forming disk. CZ08 argue that the disk only re- 
ceives the momentum of SNexp corresponding to the velocity 
vectors inclined by less than H /R. However, the fast rotating 
disk makes the SNexp anisotropic. Considering the rotation 
of the disk, the SNexp energy is channeled into two direc- 
tions: vertical and 0-direction, whose ratio is characterized by 

EL/E^n - (Vrot + Vsn)^ /Vsn =16, where Vsn « 10^\m s"' is 
the typical value of SNexp velocity (see equation 4 in Shull 
1980) and y,ot ~ 3000km s~' is the Keplerian rotation veloc- 
ity at /? = 10'*/?sch- We find E^^ > E^j^, that is to say, the 
SNexp is anisotropic in the fast rotating disk. This is the 
reason why the SNexp shell is elongated as a "banana" (Ro- 
cyzka et al. 1995). Most of SNexp energy will be exhausted 
along the 0-direction, triggering the torque to transfer angu- 
lar momentum outward. A SNexp in the SF disk is not able 
to break down the SF disk, but its evolution makes the metal- 
licity locally homogeneous. The 2D-treatment in Rozyczka et 
al. (1995) is a good approximation. The viscosity by SNexp 
is about 0.1-1. However, 3D numerical simulations should 
be carried out in the future to clarify the importance of the 
SNexp-driven viscosity in angular momentum transportation 
the disk. We use a constant a = 0.1 in the entire SF disk in 
equation (4). 

3.2. The metallicity equation 

As a natural consequence of the ongoing star formation in 
the disk as described above, the metallicity is enriching and 
evolving along with the feeding of the SMBHs. As shown by 
the cooling function in equation (25) (§4), metals play a key 
role in the cooling of the hot gas in the SF disk. In order to find 
the diffusion of hot gas, we have to calculate the metallicity 
gradients in the disk. Detailed calculations of the gradients in- 
volve many processes, including stellar evolution, metal pro- 
duction in SNexp, and the dynamics of the SF disk. Here 
we assume that stars stay at the same radius where they are 
formed until they produce SNexp in the SF disk (i.e. neglect- 
ing migration of stars as discussed in §2.1). Since the mixing 
timescale of metals with its surroundings fully depends on the 



Star Formation in Self-gravity Accretion Disks in AGNs 



TABLE 1 

A Summary of The Model of Star Forming Disk 



Parameter 


Physical meanings 


Typical values (for SMBHs with 10* M©) 


£em 


Eddington ratio 


0.3 


V 


radiative efficiency 


0.1 (for intermediately spinning SMBH) 


a 


viscosity parameter 


0.1 (from turbulence excited by the SNexp) 


MiiA 


mass of self-gravitating (or SF) disk 


~ 7rEgasi?L, ~ 5 X lO'-Mo 


SFR 


star formation rates 


~l.OM0yr-' 


SN rate 


supernova rates 


~ 3.0 X 10"^ yr"' for the Salpeter IMF 


RSG 


self-gravitating radius of the SF disk 


~ 10^i?sch ~ lO^^pc 


^in 


inner radius of the SF disk 


= Rsa 


"out 


outer radius of the SF disk 


~ sublimation radius of dusty torus ~ 1.0 pc 



SN expansion rate, we find the mixing process can rapidly ter- 
minate, in a time tj^ix "^ fsi ^ ^r^ where the radial advection 
timescale is given by Jr ^ R/Vr ^ R/aVK « 3 x lO^yr for 
a = 0.1 and R = O.lpc. The metal-rich gas then mixes with 
its surroundings, and is advected inward by the inflows while 
simultaneously being evaporated by the SNexp heating. 

A complete calculation of the metallicity evolution would 
individually follow each of the heavy elements produced 
in the SF disk, but here we follow only the overall frac- 
tion of metals relative to the progenitor mass, designated niz 
(Woosley et al. 2002). Following Wang et al. (2010), the 
fraction of the SNexp progenitor mass returned to the SF disk, 
fc = Ssn/S*, can be obtained from the initial mass function 



1 f°° 



(5) 



where M, = J m^,N^,dm^. Evidently, the factor /c depends 
on the IMF. However, the IMF in the SF disk is very poorly 
known, as even the IMF in normal galaxies is still a matter of 
debate (see the extensive reviews of the IMF by Kroupa 2007 
and Dave 2008). The state of the IMF could be a function of 
the surroundings of the star forming regions (e.g. Krumholz et 
al. 2009). However, the stellar mass function in the Galactic 
center could provide a good hint for the present case. We 
approximate the IMF as 

{NimZ for niiain < «* < »J2, 
(6) 
A^2'w* for m» > m2, 

where f3i and (32 are the indices of the IMF, m2 is the critical 
mass distinguishing the two power laws of the IMF, and the 

two normalizations are connected through N2 = Mmj for 
a continuity of the two power-laws. Top-heavy properties of 
the IMF are reflected by the two parameters /3i and m2. Usu- 
ally, (32 = 2.35 which is the Salpeter function. Figure 2 shows 
several possible IMFs and the dependence of the factors fc on 
the IMF. For a larger lower limit to the critical mass (nic) of 
stars, these approximations and assumptions work better. For 
a typical IMF, we use fc = 0.6. 

Since the lifetime of the massive stars is ~ 10''"^ yr (the 
minimum mass of stars is ^ TMq for this timescale), much 
longer than the timescale of the radial motion (fadv ^ R/Vr ^ 
R/aVK ~ lO^r for R = O.lpc in Wang et al. 2010). Massive 
stars produce SNexp in the accretion inflows, with a delay 
due to stellar evolution on a time scale of hydrogen burning 



(f,). The lag due to evolution effects must be considered for 
the time-dependent disk. Taking into account the mass used 
up in star formation, the net mass decrease due to star for- 
mation and SNexp is given by E*-/cI]',, where S, is the 
surface density of the current star formation rates at time t 
and E'^ is the rate at an earlier time offset by the lag time f,, 
namely, S',(f) = S,(f-f*). The term /cS', leads to the metal 
enrichment in the SF disk. We should note that metallicity is 
a cumulative parameter in a region with ongoing star forma- 
tion. Considering a ring AR at radius R, the gas lost to star 
formation is 27r/? A/?(E* -/^.E'^), and the advection of gas due 
to radial motion is InRARVRT^eas- We then have 



^(27r/?Sg,,) = ■^(2TrRVR^s,,)-27:R(t, 



yielding 



5S 



1 dM 



dt IttR or 



gas 



■(s,-/.s;), 



-fct'J, (7) 



(8) 



where the inflow rates M = 27r/?VREgas. We should note that 
the inflow rate varies with the radius as a result of the varying 
amount of gas used up in star formation. For a metal-enriched 
gas following equation (8), but with the inclusion of the dis- 
sipation of metals by star formation and production of metals 
by SNexp, we have 



d^z 



1 d(ZM) 



dt 2ttR or 



-Zt^ + fcniz'S'^: 



(9) 



where ZM is the metal rate taken by the inflow, and niz = 
Mz/M* is a fraction of the metal mass to its progenitor mass. 
mz G [0.1,0.4] depends on the progenitors (Woosley et al. 
2002). Introducing the metallicity Z = Sz/Egas, we have the 
metallicity equation 



dZ ainE 



dt 



■+Z- 



'gas 



1 



d(ZM) 



dt 



2ttRJ:„ 



OR 



1 

Eon 



■(ZE^-fcmzK). 

(10) 
We would like to stress that this equation can be applied to de- 
scribe the progression of AGN evolution. Provided the struc- 
ture of the SF disk is known, the metallicity and its gradient 
can be obtained from equation (10) to diagnose the evolution 
of accretion rates for both the newly-born and fading AGNs 
during a single episode. For an approximately self-similar 
evolution of accretion rates during a single episode of AGNs 
(e.g. Yu et al. 2005), such as in strong starburst galaxies com- 
posite with AGNs (e.g. Laag et al. 2010), this equation can 
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Fig. 2. — (a) Several different normalized IMFs with HI2 = 13Mq based on Kroupa (2007). (h) The parameter fc (equation 5), which is a fraction of the SNexp 
progenitor mass returned to the SF disk, for the same IMFs shown in the previous panel. For most cases, more than 60% of the gas will return to the SF disk to 
enrich the metallicity. (c) The parameter /», which is the number of SNexp per unit mass of stars formed. We use /, = 0.02Mg , a conservative value, in this 
paper. 

show the time-behavior of the metallicity. This would be an 
extremely important key to unveiling the early phase of AGN 
evolution. 

For a stationary SF disk (9Sgas/9f = 0) during a single 
episode of AGNs reaching a steady luminosity and S'^ = S*, 
we have, 

^^ ■■2TrR(l-f,)t,, (11) 



3.3. The formation of a metallicity gradient 

Radial transportation of metals due to advection produces a 
saturated metallicity gradient in the SF disk. When the metal- 
licity is saturated (i.e. dZ/dt = 0), we have 



dZ lirREt 

-—=fc(Z-mz) — -— 
dR M 



(16) 



dR 
and metallicity evolution reduces to 



from equations (12) or (13). The analytical solution is given 
by 



dZ M dZ 



■^gas 



dt 2nR dR 



+ fc(mz-Z)j:^ 



(12) 



In 



mz-Z(R) 



mz-Zo 



■■fc 



27r/?E* 



M 



-dR = fc 



Y,^dR 



(17) 



where the parameters of the SF disk are constant with time, 
except that the metallicity is increasing. A simplified version 
of equation (12) is given by 



dZ dZ E, 

^=VR^+/.(mz-Z)— , 



(13) 



which is imposed by the initial and boundary conditions 

Z(f = Q,R) = Zo, (/?i„ <R< Rout); (14) 



As an illustration of the present model for the metallicity gra- 
dient, we approximate the numerical solution (Wang et al. 
2010) as Vr = flix"', and /cE*/Egas = fl2X~"% where oi = 
2.5 X 10"^ pc yr\ fl2 = 6.5 x IQ-V"', "i = 0.94, az = 0.54, 
fc = 0.6, mz = 0.4, X = R/Roui and Rom = l.Opc (which is the 
inner radius of the dusty torus). It follows from equation (17) 



Z(x) = mz - (mz - Zq) exp 



aiRp 
aj,a\ 



■(x"^-l) 



(18) 



and 



Z(t,R = Rout) = Zo, 



(15) 



where Ri„ and /?out are the outer and inner radii of the SF 
disk, respectively. Here we assume that the initial condition 
of metallicity is homogeneous throughout the entire SF disk. 
We would like to point out that for the steady SF disk 1) 
metallicity as an integrated parameter in a steady SF disk is 
evolving, this is one of the most important characteristics of 
the present scenario; 2) the metallicity gradient and evolution 
only depend on the radial transportation if the star formation 
law is a linear relation; 3) Z{t,R) clearly relies on the param- 
eters fc and mz, but we neglect their dependence on radius 
as well as the mass of SNexp progenitors; 4) the disk struc- 
ture and metallicity actually are coupled to each other and a 
self-consistent model should produce the structure and metal- 
licity simultaneously, although we presume a known structure 
of the SF disk in this paper The largest uncertainties of the 
present results originate from the star formation law and the 
factor fc. We keep these uncertainties in mind, but leave them 
for a future paper. 



where 03 = l-ai-ax obtained from the approximation of 
disk structure. When the index term is much less than unity, 
we have the approximation 



Z(x)wZo- (mz-Zo) 



aiRo 

0:301 



■(x"^-l) 



(19) 



for the outer regions, displaying a power-law gradients as Z oc 
x"\ For a boundary with an extremely low metallicity, the 
metallicity tends to Z w mz [fl2^out(l ~ x""^) / Ci7,(i\\ ■ 

Figure 3 {left) shows the metallicity and its gradient in the 
SF disk. The gradients can be characterized by one power-law 
part as expressed by the approximation in Eq. (19) and by a 
second saturated component. The exact shape of the metal- 
licity gradient depends on the boundary condition, namely on 
the metallicity at the inner edge of the torus, but the inner- 
most region of the SF disk always has the maximum metal- 
licity independent of this boundary condition. The maximum 
metallicity is Z^ax = niz, which is the metallicity of the SNexp 
progenitor. 



Star Formation in Self-gravity Accretion Disks in AGNs 



10 



'3 



-1-1 



^ 

N 



0.1 




z/z 



0.01 




0.01 



Time (10^ yr) 



0.1 1 



0.1 
Radius (1 pc) 



Fig. 3. — Left Metallicity gradient of the SF disk wlien Z is saturated at each radius. The four colored Hnes are for different initial metallicity in the SF disk. 
We find that the gradient of the metallicity is quite large, so there should be significant difference between the metallicities of the inner and outer parts of the 
accretion disk. Right Metallicity evolution of the SF disk for a given initial metallicity. The model of metal evolution is only valid after a few times lO^yr after 
the initial star formation. We use M, = IO^Mq and ^ej,] = 0.3. 



3.4. Metallicity evolution 

If a SMBH is fed with a constant accretion rate in a single 
episode, its evolutionary progress cannot be probed by using 
the SMBH mass accretion rate, but the metallicity is an indi- 
cator of the ongoing star formation. This interesting feature 
allows us to estimate the age of the SMBH activity. The par- 
tial differential equation (13) can be analytically solved for 
the case of a SF disk with a power-law structure. Employ- 
ing the approximation of Wang et al. (2010), Vr = oix"', 
/cS*/Sgas = a2X~°'-, the solution to the partial differential 
equation (13) is characterized as 



+ — !-f = Ci, 

l-Qfl Rout 



and 



In(mz-Z) 



fl2^o 



fllQ;3 



= C2, 



(20) 



(21) 



where a^ = 1 —a\ — a2, C\ and C^ are two constants. We thus 
have the general solution for the metal evolution as 



ln{mz-Z) 



aia^ 



= T 



+ 1 

\-a\ Rout 



(22) 



where T is any continuous function. Applying the general 
solution to the initial and outer boundary conditions, we have 
the solutions 



In 



aiRo 



/ mz-Z 
\mz-ZoJ aia3 



"out 



and 



In 



mz-Z 

mz-Z) 



aiRo 



a\a-i 



-(jc"-'-l). 



aj/Cl-Ql)^ 

(23) ' 
(24) 



For a given SF disk with typical parameters. Figure 3 right 
shows the evolution of the metallicity gradient. At each ra- 
dius in the SF disk, the metallicity increases with time until it 
reaches a saturation value at a critical time. These saturated 



points consist of a breaking curve on the metallicity surface. 
Radial transportation due to advection balances with the lo- 
cal production of metals. In the innermost regions of the SF 
disk, the metallicity reaches its maximum value in a few times 
10^ yr It should be noted that the radial transportation is due 
to the SNexp-turbulent viscosity (Wang et al. 2010). The 
metallicity gradient itself is thus a self-organized result with- 
out external factors. The property of metallicity with time is 
very useful for estimation of star forming disk ages. 

Equation 10 describes the metallicity evolution in the disk, 
provided that the mass rates of the supplied gas are known. It 
is clear that galaxies are undergo repeated episodes of SMBH 
activity and star formation (e.g. Marconi et al. 2004; Wang 
et al. 2006; also see the introduction in CZ08). For a steady 
star forming disk, the equation reduces to one for metallic- 
ity evolution in a single episode of SMBH activity. During 
that single episode, the evolutionary phase of SMBH accre- 
tion is very hard to estimate from just the accretion rates since 
they are roughly constant. However, the property that metal- 
licity monotonically increases with time allows us to use it 
as a probe of the age of the AGN outburst. The next sec- 
tion will discuss for a single episode the diffusion of hot gas 
driven by SNexp from a steady star forming disk, forming a 
hot "corona" above the disk. 

4. THERMAL DIFFUSION OF HOT GAS FROM THE SF DISK 

Interaction between the blast wave of a supernova rem- 
nant and interstellar clouds will lead to the evaporation of 
clumps, producing hot gas between the clumps (McKee & 
Cowie 1975), meaning that SNexp heating is driving produc- 
tion of hot gas in the SF disk. A stationary evaporation rate 
from the SF disk can be reached by the balance between the 
SNexp heating and the cooling of the hot gas. Metallicity 
has a strong influence on the cooling rate of hot plasma (e.g. 
Bohringer & Hensler 1989). Since metallicity is quite high 
in AGNs and quasars, we include the effects of metallicity on 
the cooling of the hot plasma, and get the diffusion of hot gas 
self-consistently. 

With the inclusion of emission lines, the cooUng rate due to 



WANG ET AL. 



thermal emission from hot gas is approximately given by 



A.^^„, = 2x10- 



■10 27,-1.3 



erg s ' cm ^, 



(25) 



for a range of temperature of 10^ - lO^K from Fig.l in 
Bohringer & Hensler (1989). Here Te = T/IO^K and tif, = 
tie/lQi^ cm"^ are the temperature and density of the hot gas 
evaporated in the SF disk. 

The heating rate for the same gas depends on the rate of 
SNexp and the energy injected per SNexp. Using the IMF 
given by equation (6), the number of SNexp per unit mass of 
stars formed is 

1 f^ 
f* = T7- N,dm,. (26) 

Values of /* depends on the IMF, but it is in a range of 
/* G [0.02,0.05]Mq from Figure 2c. We use the lower limit 
value of/* = 0.02Mq in this paper. The heating rate per unit 
volume due to SNexp is then given by 



ASN =/*(;^l^'^sN 



= 1.0 X 10-iVo.02ei/6S*^o//r4'^5i erg s"' cm'^, 



(27) 



where /0.02 = /*/0.02Mq, ^i/e = C/| is the fraction of 
the kinetic energy converted into thermal energy, E* = 
S*/10"MQyr"'pc~^, £51 =£'sN/10^'erg is the kinetic energy 
of the SNexp and H14 = i/a/lO'^cm. The values of ^ and 
/, are both uncertain, and are interrelated. In the following 
equations they appear together as the product /^ = ^/, . For 
a given /^, the higher /,, the smaller ^. For the minimum 
/*.min = 0.02Mq , the parameter is ^ « 1/6 roughly. Provided 
/^ is at a level of ~ 1 0~^Mq , which is sufficient for the present 
model, the required ^ is only a few percent. Though it is diffi- 
cult to estimate /» and ^, the present model is not sensitive to 
their exact values. 

We can then estimate the temperature of the hot gas from 
the energy balance A^^^j = Asn in the SF disk. This yields 



nJT.-'-^Zi = 0.50 /4,-3£5iS,,oHri, 



(28) 



where Zi = Z/Zq and /^__3 = /^/3 x lO"-'. A further constraint 
on the condition of the hot gas is that it should not be cooled 
before escaping from the SF disk. This requires a critical 
value of f^j^, = Ha/cs, where Cj = {kT / iinip)^/^ is the sound 
speed of the hot gas, and 



'cool 



riAT 1 



(7- Dm A; 



SP =3.1xloVo'«6'762X'sec, (29) 

cool 



where k is the Boltzmann constant, nip is the proton mass, 
7 = 4/3. We then have 



3/2, 



nl'Tl^Z\' =2.91 iil^^H,^ 



(30) 



Combining equations (28) and (30), we have the density and 
temperature of the hot plasma in the SF disk, 

„ _n 00 /-O.eSpO.feS, 0.45 7-0.36f,0.65r,-0.36 /oin 

n6-0.88/j_3/i51 ^0 ^1 ^*,0 ^14 : (3i) 



and 



Since the hot gas is cooling on a time scale longer than the 
timescale for escaping from the SF disk, it forms strong diffu- 
sion into regions above the SF disk. The diffusion rates driven 
by the SNexp are then given by 



Miifi = / 2'!TCsmpniiRdR 

'Kin 



1.8 / CoZf-^*t%'H^-^'x'dx'MQyr-\ 



(33) 



where Co = /IZg^si^A^o ^'. R is the radius of the SF disk, x' = 
R/lO~°-^pc, Rin = RsG and/?out= Ipc. Equations (31), (32) and 
(33) are the form solution of the hot gas in the SF disk in light 
of S,, //d and the metallicity Z. 

Inserting equations (3), (4), (18) and (22) into (31-33), we 
find the structures of the diffused hot gas for the case of sat- 
urated metallicity and the evolution of those structures with 
increasing metallicity, respectively. For the saturated metal- 
licity. Figure 4 left shows the structure of the diffused hot 
gas above the SF disk. The left panel shows that the den- 
sity and temperature distribution of the hot gas is insensitive 
to the boundary metallicity when the metallicity reaches its 
saturated state. The temperature of the hot gas in the SF disk 
is roughly constant with radius, staying in the range between 
(0.5 ^ 1 .0) X lO^K. We note that the temperature is lower than 
the typical Compton temperature of the AGN gas. The den- 
sity of the hot gas strongly depends on the radius, falling from 
roughly 10^ to 10"*cm"^ from the inner to the outer part of the 
disk. We find that the total mass rates of the hot gas diffu- 
sion over the entire SF disk can be as high as a few Moyr"' 
if the Eddington ratio is i?Edd = 0.3. We should note that this 
mass rate could be the lower limit since we have conserva- 
tively taken /^ ~ 10~^Mq . Additionally the diffusion rates 
of hot gas are comparable with or exceed the Eddington limit 
of a IO^Mq SMBH, but they are still smaller than the rates 
of the inflows developed from the inner edge of a dusty torus 
(see Wang et al. 2010). 

Figure 4 right panel shows the evolution of the hot gas and 
its diffusion with the metallicity. We find that temperature 
and the density of the hot gas vary only slightly with Z, but 
the diffusion rates vary by a factor of a few. The higher the 
metallicity, the lower the diffusion rates. For an individual 
AGN evolving with metallicity, the thermal diffusion becomes 
weaker. 

It should be pointed out that the diffused hot gas is actu- 
ally still bounded by the potential of the SMBH. We can esti- 
mate the thickness of the diffused hot gas by assuming verti- 
cal static balance, //diff « cjn^ = 1.4 x 10^'' T^^^rl^^Mgcm, 

Hmf/R ^ 0.05 Tg / Vy^ and //diff « 3//d, where Ha is the thick- 
ness of the SF disk. The thermal diffused hot gas then meets 
complicated fates as a result of heating by radiation from the 
accretion disk. The hot ionized gas will undergo thermal in- 
stability and form cold clouds in the broad line region since 
it is continuously supplied through the evaporation driven by 
SNexp in the SF disk. The metallicity gradient in the SF disk 
determines the metallicity in the BLR. A detailed model of 
the BLR formed by the hot gas will be the subject of a forth- 
coming paper (Wang et al. 201 1 in preparation). 
Combining equations (3), (4) and (33), we find that the dif- 
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Fig. 4. — Structures of the hot gas diffused from the SF disk. Left is for an steady SF disk and right for the time-dependent. From bottom to top panels, we 
show temperature, density, surface-density of star formation rate and mass injection rate of the thermal diffusion. We use the Eddington ratio £em = 0.3. The 
mass injection rate is from the [i?i„,S] part of the SF disk, the [i?in,Rout] (i.e. Equation 33) is the total mass rates with an order lM0yr"' . We assume the initial 
metallicity Z = Q.SZq throughout the SF disk (the rgiht panel). 



important observational implications, indicating that the dif- 
fusion of hot gas is more sensitive to the SMBH mass than 
metallicity and Eddington ratios. In light of Z ex L oc M.^Edd 
from observations (Hamann & Ferland 1999), we have Mdiff oc 
^Edd^»'^^- ^^ is therefore expected that the diffusion of the 
hot gas is not sensitive to the Eddington ratio of the accreting 
SMBH, but is moderately sensitive to the SMBH mass. For 
an individual quasar, Minj oc Z""^"*, implies that quasars in the 
late phase of a single episode should have relatively weaker 
outflows, but higher metallicity. Last, the diffusion rates of 
the hot gas may be greatly enhanced if /^ = ^/^ > IQ~^Mq 
for a more top-heavy SF disk. It should be kept in mind that 
equation (32) is a conservative estimate of the diffused hot gas 
from the disk. 

We note that the hot gas inside the SF disk may have soft X- 
ray emission. Total emission from the hot gas can be simply 
estimated by the free-free emission in the disk. We have Luc = 
J 2nj(ne, Te)RHdR « 2.0 x 10"*" / r^'^-^^dr^ « 5.0 x lO""' erg/s, 
where the emissivity of the free-free emission is y « 1 .4 x 

lO'^^HgTe erg s"' cm"-'. The free-free emission in soft X- 
rays and radio would be negligible. It should be noted that 



in our model the thermal diffusion is driven by a different 
mechanism than in the model of winds driven by starbursts 
(Chevalier & Clegg 1985). The present model stresses the role 
of evaporation of the molecular clouds through their heating. 
In the starburst model, bubbles are energized as winds by the 
kinetic energy of the starburst in a physically much thinner 
medium than the SF disk. 

Finally the influence of the diffusion of hot gas on the disk 
can be estimated from the ratio of the disk accretion rate to 
the mass rate of the diffusion. We find that Mdiff ~ IMoyr"' 
is significantly smaller than the inflow mass rates (Mi„flow ^ 
lOM0yr"') (see Figure 2 in Wang et al. 2010). The influence 
of the diffusion of hot gas on the disk can be generally ne- 
glected in the star forming regions, but it might be important 
in the region around the self-gravity radius, where the diffu- 
sion mass rates are comparable to the accretion rates onto the 
SMBH. In this region, the complicated interaction between 
the warm gas and the accretion flow onto SMBHs may lead 
to some observable features in X-ray spectra. We leave this 
for a future investigation. On the other hand, cooling of the 
hot winds may be balanced by the efficient Compton heat- 
ing driven by the X-ray photons emitted from the hot corona 
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above a Shakura-Sunyaev accretion disk. The final fate of 
the outflowing hot gas is a very complicated question, since 
this gas will be subject to a complex array of processes as it 
moves out into the host galaxy. This will be the main topic 
of a follow-up paper, but in the next section we outline some 
observable effects on the BLR gas. 

5. OBSERVATIONAL CONSEQUENCES FOR THE BLR 

The BLR gas is assumed to be supplied by diffusion of the 
hot gas. There can be a factor of ten or more metallicity dif- 
ference between the inner and the outer parts of the SF disk, 
leading to observable consequences. This gas is then pho- 
toionized by the radiation from the accretion disk. We can 
expect that the outer regions of the BLR are fed by the outer 
part of the SF disk. Since this gas is farthest from the contin- 
uum source, it will have a relatively low value of the ioniza- 
tion parameter, defined by L'^ = Lio„/4'!TR^c{e)n, where Lion is 
the ionizing luminosity, (e) is the averaged energy of ionizing 
photons and n is the density and temperature of the ionized 
gas, respectively. Emission lines with low ionization energy 
will be produced from this low ionization line (LIL) region. 

The inner regions of the SF disk, being closer to the ionizing 
source, will have a higher value of U, and thus produce high 
ionization lines (HILs). It should thus be possible to exam- 
ine the metallicity gradients through comparison of the LlLs 
and HILs. Additionally, measuring metallicity differences be- 
tween the HIL and LIL regions can provide a rough estimate 
of the age of the single episode of the SMBH activity, and po- 
tentially even reveal its evolutionary track. It is thus expected 
that the metallicity measured from the LlLs should be lower 
than the HILs. In principle, this theoretical prediction can be 
examined by separate measurements of the metallicity in the 
HIL and LIL regions. 

Generally, metallicity in AGNs is much higher than in their 
host galaxies. This is true even after extrapolating the ob- 
served metallicity gradient in galaxies to the innermost re- 
gions. However, it is much more difficult to measure the 
metallicity in AGNs than in stars. The original way to mea- 
sure AGN metallicity (Hamann & Ferland 1993) was to use 
HIL ratios of N v/C IV or He II /C IV together with de- 
tailed photoionization models. Later, broad-LlL ratios such 
as, N III]/C III] and N lll]/0 III] were used together with 
the HIL lines to estimate the metallicity based on the LOG 
model (Hamann et al. 2000; Warner et al. 2003, but see the 
original suggestion of Shields 1976). Warner et al. (2003, 
see Figure 10; and 2004) found that the metallicity measured 
by N III]/C III] is systematically lower than N v/C IV by 
a factor of a few for a sample of ~ 800 AGNs and quasars. 
Although they did not realize that these differences could be 
intrinsic and instead took the averaged values to be the over- 
all BLR metallicity, their measurements agree well with the 
factor ^10 metallicity gradient across the SF disk that we 
find here. Assuming that both kinds of metallicity indicators 
properly reflect the real metallicity, the systematic difference 
found by Warner et al. (2003; 2004) clearly supports the idea 
that there is gradient of star formation within the central Ipc 
region. 

Additionally, the detailed nature of the broad line regions 
is still very poorly understood. However, the present simple 
model can conveniently explain the significant difference of 
metallicity represented by N v/C IV and N Ill]/C III] (Warner 
et al. 2003; 2004). Further high quality spectra of quasars 
are needed to measure the metallicity accurately enough to 
determine the quasar's evolutionary status. 



We would like to point out that there should be little or 
no mixing of metals between the LIL and HIL regions, even 
though the dynamical time scales in the HIL and LIL are only 
a few tens of years. In a new model of the BLR based on the 
present SF disk, we will show (Wang et al. 2011 in prepara- 
tion) that there is a mass circulation between the BLR and the 
SF disk, keeping a constant mass in the LIL regions, whereas 
the cold clouds in the HIL could be blown away by the radi- 
ation pressure. The two regions maintain a stationary state, 
showing different metallicity in LIL and HIL regions. 

6. CONCLUSIONS AND DISCUSSIONS 

In this paper, we have extended the study by Wang et al. 
(2010) of the properties of star forming disks within the cen- 
tral 1 pc regions of quasars and AGNs, where star forma- 
tion can drive an adequate fuel supply to the SMBH through 
SNexp. We have made progress in the following areas: 

• A metallicity equation is derived, unveiling a metallic- 
ity gradient in the star forming disk. We find that metal- 
licity can reach up to ^ 1O-2OZ0 in the innermost re- 
gions of the disk and a few times solar abundance in the 
outer part. As an integrated physical parameter, metal- 
licity is evolving even in disks that dynamically are in 
a steady state, and the metallicity can reach saturation 
in a few lO^years. Metallicity can be used as a "clock" 
for indicating the age of a star forming disk during one 
episode of SMBH activity. 

• Thermal diffusion with a mass rate of ^ IMoyr"' in- 
evitably developes from the star forming regions. This 
rate is comparable to the outflow rates estimated from 
X-ray spectroscopic observations. The hot gas diffused 
above the disk is likely to be the source of the observed 
broad-line-region gas. 

• The BLR should then share the metallicity gradient of 
the SF disk. The metallicity gradient might already 
have been measured by the systematic difference in the 
metallicity measured by the low (e.g. N Ill]/C III]) and 
high (N v/C IV) ionzation lines in AGNs and quasars. 

The properties of the star forming disk have important the- 
oretical and observational consequences. It is well known 
that AGNs have many separate ingredients, including the BLR 
(with separate HIL and LIL regions), the NLR, the accretion 
disk, the self-gravitating disk and the torus. However, the ori- 
gins of these different components as well as the physical con- 
nections among them are extremely poorly understood. Since 
the diffusion rates of the hot gas evaporated from the SF disk 
are comparable with the accretion rates of the SMBHs, the 
presence of this gas will have significant consequences. The 
ionized gas exposed to the AGN irradiation will have a very 
complicated fate in light of AGN heating, cooling and cou- 
pling with its dynamics. However, we speculate that the ac- 
cumulation of the diffused gas continually supplied from the 
SF disk can drive the formation of BLRs through dynamical 
and thermal instabilities. Dependence of the properties of the 
diffused hot gas on the radius would then lead to onion-like 
structures within the broad line regions. 

It is often suggested that broad line regions, broad absorp- 
tion lines, and "warmer" absorbers are the result of outflows 
(e.g. Murray & Chiang 1997). The present scenario pro- 
vides a gas supply to the broad line regions, warmer absorbers 
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etc. This points towards a unified model in which broad line 
clouds, outflowing clouds and warmer absorbers are intrinsi- 
cally linked with the star forming disk, and hence with the 
accretion disk around the SMBHs. Such a model will be pre- 
sented in a forthcoming paper (Wang et al. 201 1 in prepara- 
tion). 

A direct proof of the presence of massive stars in the core 
of a SF disk would be the detection of the spectral features 
of red supergiants. In the near future, ALMA (Atacama 
Large Millimeter Array) spectra of galactic nuclei should 
have the sensitivity to detect the CO bandhead in the near IR 
that is predicted by our scenario. Another approach would 
be to use powerful instruments to search for Ca II triplet 
(AA8498, 8542, 8662A) absorption from the nuclear regions. 

Finally, newly born AGNs and quasars may have rapidly 
growing accretion rates with self-similar behaviors such as 
Sgas oc R^'f- (Yu et al. 2005, but see Schawinski et al. 2010). 
This would allow us to obtain the general solution for the 
metallicity evolution during the early stage of the AGN phase. 



Interestingly, the metal-poor quasars have metallicities sim- 
ilar to their host galaxies (Z < O.6Z0), but they are only a 
small fraction of quasars (Fu & Stockton 2007). As has been 
suggested by Silk & Rees (1998) and others, Ultra-Luminous 
Infrared Galaxies (ULIRGs) with higher star formation rates 
(e.g. Laag et al. 2010) and low metallicity are likely to be 
in the early stage of AGN activity, soon after their birth. Our 
work here helps move us into a position to apply much more 
quantitative tests to this general idea. 

We are grateful to the anonymous referee for a very useful 
report improving the manuscript. H. Netzer and A. Laor are 
thanked for helpful suggestions and comments. We appreci- 
ate the stimulating discussions among the members of IHEP 
AGN group. Conversations with J.-C. Wang are appreciated. 
The research is supported by NSFC-10733010and 10821061. 
JAB is grateful to J.-M. Wang and the Institute of High En- 
ergy Physics of the Chinese Academy of Sciences for their 
very gracious hospitaUty, and to NASA grant NNX10AD05G 
for financial support. 



REFERENCES 



Artymowicz, P., Lin, D. N. C. & Wampler, E. J. 1993, ApJ, 409, 592 

Bohringer, H. & Hensler, G. 1989, A&A, 215, 147 

Brocksopp, C, et al. 2006, MNRAS, 366, 953 

Brunner, H., et al. 1997, A&A, 326, 885 

Chen, Y.-M., Wang, J.-M., et al. 2009, ApJ, 695, L130 

Chevalier, R. A. & Clegg, A. W. 1985, Nature, 317, 44 

Collin, S. 1999, Phy. Rep., 311, 463 

Collin, S. & Hure, J.-M. 1999, A&A, 341, 385 

Collin, S. & Zahn, J. P. 1999, A&A, 344, 433 (CZ99) 

Collin, S. & Zahn, J. R 2008, A&A, 477, 419 (CZ08) 

Dors, O. L., Storchi-Bergmann, T., Riffel, R. A. & Schmidt, A. A., 2008, 

A&A, 482, 59 
Dave, R. 2008, MNRAS, 385, 147 
Duschl, W. J. & Britsch, M. 2006, ApJ, 653, L89 
Fu, H. & Stockton, A. 2007, ApJ, 664, L75 
Gammie, R. 2001, ApJ, 553, 174 
Goodman, J. 2003, MNRAS, 339, 937 
Hamann, F W. & Ferland, G. J. 1993, ApJ, 418, 1 1 
Hamann, F W. & Ferland, G. J. 1999, ARA&A, 37, 487 
Hamann, F W., Netzer, H. & Shields, J. C. 2000, ApJ, 536, 101 
Hobbs, A., Nayakshin, S., Power, C. & King, A. 2011, MNRAS, 18333.x 
Hure, J.-M. et al. 1994, A&A, 290, 34 
Jiang, Y.-F & Goodman, J. 201 1, ApJ, 730, 45 
Kennicutt, R. C, Jr. 1998, ARA&A, 36, 189 
Kennicutt, R. C, Jr. et al. 2007, ApJ, 671, 333 
Kolykhalov, P I. & Sunyaev, R. A. 1980, Sov. Astr Lett., 6, 357 
Kroupa, R, 2007, in Valls-Gabaud, D., Chavez, M., eds, ASP Conf. Ser., 

Resolved Stellar Populations, in press (astro-ph/0703124 1 
Krumholz, M. R., McKee, C. F & Tumlinson, J. 2009, ApJ, 693, 21 6 
Laag, E., Croft, S. Canalizo, G. & Lacy, M. 2010, arXivrrOTO.17041 
Laor, A. & Netzer, H. 1989, MNRAS, 238, 897 
Lin, D. N. C. & Pringle, J. 1987, MNRAS, 225, 607 
Lu, X.-R. 2008, ChJA&A, 8, 50 
Lynden-Bell, D. 1969, Nature, 223, 690 
Malkan, M. A. 1983, ApJ, 268, 582 
Marconi, A., et al. 2004, MNRAS, 351, 1390 
Matsuoka, K., Nagao, T., Marconi, A., Maiolino, R. & Taniguchi, Y., 2011, 

A&A, 527, 100 



McKee, C. F & Cowie, L. L. 1975, ApJ, 195, 715 
McKee, C. F & Ostriker, E. C. 2007, ARA&A, 45, 565 
Murray, N. & Chiang, J. 1997, ApJ, 474, 91 
Nayakshin, S. & Sunyaev, R. 2005, MNRAS, 364, 223 
Nayakshin, S. & Hobbs, A. 2008, JPhCS, 131, a2011 
Onodera, S., et al. 2010, arXiv 1009. 1971 
Paczynski, B. 1978, Acta Astr. 28, 91 

Rafokiv, R. 2009, ApJ, 704, 28 1 

Rahman, N., et al. 2010, arXiv: 1009.3272' 

Rozyczka, M., Bodenheimer, P & Lin, D. N. C. 1995, MNRAS, 276, 1179 

Sakimoto, P J. & Coroniti, F V. 1981, ApJ, 247, 19 

Schawinski, K., et al. 2010, arXiv: 101 1.0427 

Shakura, N. Z. & Sunyaev, R. A. 1973, A&A, 24, 337 

Shields, G. A. 1976, ApJ, 204, 330 

Shields, G. A. 1978, Nature, 272, 706 

Shlosman, I. & Begelman, M. C. 1987, Nature, 329, 810 

Shlosman, I., et al. 1989, Nature, 338, 45 

Shlosman, I. & Begelamn, M. C. 1989, ApJ, 341, 685 

Shull, J. M. 1980, ApJ, 237, 769 

Silk, J. & Rees, M.J. 1998, A&A, 331, LI 

Simon, L. E. & Hamann, F 2010, MNRAS, 407, 1826 

Smith, M. & Mac Low, M.-M. 1997, A&A, 326, 801 

Sun, W. H. & Malkan, M. A. 1989, ApJ, 346, 983 

Tan, J. & Blackmann, E. 2005, MNRAS, 362, 983 

Thompson, T., Quataert, E. & Murray, N. 2005, ApJ, 630, 167 

Wang, J.-M., Chen, Y.-M. & Zhang, F 2006, ApJ, 647, L17 

Wang, J.-M., Yan, C.-S., Gao, H.-Q., Hu, C, Li, Y.-R. & Zhang, S. 2010, 

ApJ, 719, L148 
Warner, C, et al. 2003, ApJ, 596, 72 

Warner, C, Hanmann, F & Dietrich, M. 2004, ApJ, 608, 136 
Wheeler, J. C, Mazurek, T. J. & Sivaramakrishnan, A. 1980, ApJ, 237, 781 
Woosley, S. E., Heger, A. & Weaver, T. A. 2002, Rev. Modem Phys., 74, 

1015 
Yu, Q., Lu, Y.-J. & Kauffmann, G. 2005, ApJ, 634, 901 



